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I Abstract 

When used with coherent hght, optical imaging systems, even diffraction-hmited, are inherently unable to reproduce both the 
C/3 (implitude and the phase of a two-dimensional field distribution because their impulse response function varies slowly from point 
^ to point (a property known as non-isoplanatism) . For sufHciently small objects, this usually results in a phase distortion and has 
-*-^ ho impact on the measured intensity. Here, we show that the intensity distribution can also be dramatically distorted when objects 
^~bf large extension or of special shapes are imaged. We illustrate the problem using two simple examples: the pinhole camera and 
, the aberration-free thin lens. The effects predicted by our theorical analysis are also confirmed by experimental observations. 
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'^PACS: 42.30.Kq, 42.30. Va, 42.30.Lr 
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1. Introduction 
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Current technology, especially the ability to manufac- 
ture aspherical surfaces, allows lenses and mirrors to be de- 
signed that minimize the most important geometrical aber- 
rations. Such optical elements are nearly ideal instruments 
obeying the laws of Gaussian optics even for far off-axis 
points and non-paraxial rays. Self-luminous objects or ob- 
jects illuminated with incoherent light can be imaged with 
high fidelity. The quality of the resulting image is merely 
fixed by the resolution of the instrument which is related to 
its numerical aperture. The instrument itself is said to be 
'diffraction-limited and can be considered as a linear filter 
for the intensity of light |1|2| . 

When imaging objects with coherent light, the condi- 
tions for accurate image formation are more severe since 
both the relative amplitudes and the relative phases of the 
object points have to be mapped to the corresponding im- 
age points (up to the resolution capability of the instru- 
ment). This only happens if the response of the optical in- 
strument to the field from a given point source is indepen- 
dent of its position in the object plane, or in other words, if 
the coherent impulse response |1I2J of the system is space- 
independent. The instrument then acts as a linear filter for 
the complex field amplitude. According to the terminology 
of [1] , such an instrument is said to be isoplanatic. In gen- 
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eral, even aberration-free optical instruments designed to 
map some planar object to an image plane under incoher- 
ent illumination do not meet this condition. Some spatial 
phase distortion is unavoidably introduced, which, when 
combined with a finite resolution, severely modifies the in- 
tensity distribution of the image. This was first recognized 
by Dumontet [7 , and later by Tichenor and Goodman [1] , 
who showed that a thin lens can only be considered as an 
isoplanatic system if both the object and the image lie on 
spherical surfaces So and Si, tangent to the geometrical- 
optics object and image planes O and X respectively, and 
having their center of curvature in the plane of the lens (see 
SecO. 

In practice, an aberration-free optical instrument can be 
treated as an ideal coherent-light imaging system whenever 
the spherical surfaces So and Si can be approximated by 
their tangent planes O and X. We emphasis that this is only 
viable when the object to be imaged is very small and lies 
close to the optical axis. A weaker imaging condition has 
been obtained by Tichenor and Goodman who showed that 
non-isoplanatism of a thin lens has a negligeable effect on 
the intensity distribution of the image if the object diameter 
is smaller than about a quarter of the lens diameter [?j . In 
this paper, we investigate the effect of non-isoplanatism on 
coherent image formation when this condition is not met. 
That situation may be encountered in many fields of op- 
tics where large-sized objects are imaged through powerful 
limited-aperture instruments, as in coherent far-field mi- 
croscopy, optical lithography [5] , holographic data-storage 
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[5] , and dipole-trapping of neutral atoms \1\ . In particular 
arrays of coherently emitting point sources, like individual 
trapped atoms excited by the same laser beam [8191, will 
be subject to the phenomena discussed here. Effects simi- 
lar to those that we described in the context of imaging are 
also expected when lenses are used to perform the spatial 
Fourier transform of a two-dimensionnal field distribution, 
as or holographic dipole-trapping of atoms [10] . 

2. Non-isoplanatism of a pinhole camera 

The problem of non-isoplanatism can be best understood 
by considering the simple example of coherent image for- 
mation through a pinhole camera (Fig. [T]). 

Consider the following situation: Two mutually coherent 
point-objects, Pi and P2, are imaged through a small pin- 
hole (of radius a) from an object plane O to an image plane 
X. We call {xo.k,yo.k) the transverse coordinates of the point 
sources Pk {k £ {1, 2}), and {xi^k,yi,k) the transverse coor- 
dinates of their geometrical images P^ . For simplicity, we 
chose yo.k = Vi.k = and assume that the point sources 
emit in phase. We assume that the pinhole is so small that 
the resolution is limited by diffraction. The "images" of Pi 
and P2 are two Airy patterns centered on the geometrical 
image points P{ and P2- Let's consider that P{ and P2 are 
at the resolution limit according to the Rayleigh criterion 
Xi^i — Xi^2 — 0.61Azi/a. Due to the geometry displayed in 
Fig. [U the relative phase of their respective Airy patterns 
is 4> = 27r(ri — r2)/A. This shows that the phase relation 
between points is not mapped properly from the object plane 
to the image plane. Furthermore, the relative phase (p varies 
as a function of Xo^m, = ixo,i + Xo,2)/'2, the mean distance 
of the point sources to the optical axis. The phase (j) is obvi- 
ously null if Xm — 0, but it already reaches the value tt/2 
when Ixo.ml — a/2.44. Since a is usually about 1 mm for 
a pinhole camera, (j) varies extremely rapidly when the two 
point sources are moved over the object plane. Fig.[2]shows 
the intensity distribution in the image plane resulting from 
the interference of two Airy patterns with relative phase 
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Fig. 1. Scheme for coherent image formation by a pinhole camera. 
The points P^ and P^ in the image plane X are the geometrical images 
of the point objects Pi and P2 in the object plane O. The light-gray 
zones around P[ and Pj ™ the X-plane represent the individual 
diffraction patterns resulting from illumination of the pinhole by Pi 
and P2 respectively. The dark-gray zone represents the area where 
the individual diffraction patterns overlap and interfere. 
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Fig. 2. Intensity profile in the image plane of a diffraction-limited 
pinhole camera due to the interference of the Airy patterns of two co- 
herent point sources. The distance \xo,i — Xo.2\ between the sources 
is equal to (a) 0.61Xzo/a (the Rayleigh resolution limit), and (b) 
0.915Xzo/a (1.5 times the Rayleigh resolution limit). The Airy pat- 
terns have been normalized so that their peak intensity is 1. Panel 
(c) and (d) show the intensity profiles along the a;-axis corresponding 
to the images in panel (a) and (b), respectively. Vertical arrows in 
panels (c) and (d) point the positions of the centers of the interfering 
Airy patterns. 

(f> — 0, 7r/2, and tt. In panel (a) the distance between the 
sources corresponds to the Rayleigh resolution criterion, as 
discussed before; in panel (b) this distance is increased by 
a factor of 1.5. Panel (c) and (d) display the intensity pro- 
file along the x-axis for images in panel (a) and (b) , respec- 
tively. Note that the intensity distribution for </> = tt/2 is 
the same as for incoherent sources. Because of this inter- 
ference, the image points may be unresolved even if their 
relative positions fulfill the Rayleigh criterion. 

The preceding example shows that different intensity dis- 
tributions must be expected from identical object patterns 
depending on their positions in the object plane. This is 
due to an incorrect phase mapping between the object and 
image plane. The relative phase of an object point and 
its image depends on the position of the point source in 
the object plane. We refer to this situation as the non- 
isoplanatism of coherent imaging. Considering the pinhole 
camera as a linear optical system [5], non-isoplanatism is 
related to the space- variance of its impulse response. We 
develop this point of view hereafter. 

Since the pinhole camera is a linear system, the complex 
field amplitudes in the object and image planes, Uo{xo, yo) 
and Ui{xi,yi) satisfy the integral relation 
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{xi,yi) = JJ h{xt,yt\xo,yo) Uo{xo,yo)dxodyo. 



The impulse response of the camera (in the paraxial and 
far-field approximation) is given by the Fraunhofcr integral 
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where {xo, yo) are the coordinates of the point source in the 
object plane, (x^, yi) are the coordinates of the "observation 
point" in the image plane. The integration domain A is the 
pinhole opening disc, while 

?' = ^/xfTyfT^f , 

^1 + yf + zf 

are the distances from the point source and the observation 
point to the pinhole center, and 

1 /a;. 
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are the spatial frequencies of the plane waves diffracted by 
the pinhole. On the second line of ([T]), the function 
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has been introduced {Ji{x) is a Bessel function of first 
kind). Its square S'ff{x, y) is the usual intensity point spread 
function or Airy pattern of diffraction-limited optical sys- 
tems and it is normalized so that \\n\\i^\^_,^5isi{x,y) = 
S{x,y) (the Dirac distribution). The "scaling" property, 
5N{Cx,Cy) = ScN{x,y)/C^ for any C G Rq, can be used 



to simpfify ([T]): 

h{x,,y^\xo,yo) = \M\ e^*('^+^) x 



6^{xi - Mxo, yt - Myo), 



(3) 



where M = —sjr is the geometric magnification ratio of the 
camera. Note that we can make the approximation M « 
—Zi/zo because r and s vary only very slightly over the ob- 
ject and image planes. For the same reason, ^^/(As) {x, y) ~ 

Sa/{\z-){x,y)- 

The phase factor in Eq. ([3]) plays an important role. 
If the phase factor were not there, the impulse function 
corresponding to any point {Xq, yo) lying in the object 
plane would be the same as the impulse function of the 
origin (0, 0), but translated to the geometrical image point 
{Mxo,Myo); in that case, the system would be space- 
invariant, or isoplanatic, and the image of a coherent ob- 
ject with field amplitude Uo{xo, yo) could be computed by 
simply convoluting Uo{xo,yo) with the impulse function. 
In other words, the pinhole camera would act as a linear 



filter that reproduces the objet Uo{xo,yo) with a lower 
resolution. The example of Fig. [2] clearly shows that this 
is not the case. Because of the phase factor, the impulse 
response function is different for different points in the ob- 
ject plane. This creates interference patterns that modify 
the image more significantly than a simple blur. The image 
of a field profile Uo{xo, yo) is given by 
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To establish Eq. ^ , we make use of the fact that 
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(5_a. {x^ - Mxo, yi - Myo) = j^5^ (^ " ^°' M ^ ^°^ 

can be considered either as a function of (x^, yi), the usual 
Airy pattern centered on the geometrical image point 
{Mxo, Myo), or as a function of {xo,yo) centered on the 
object point {xi/M,yi/M). Note that in the first inter- 
pretation, the first-zero full-width of the Airy pattern 
Ci = 1.22Xzi/a represents the region of the image plane 
most influenced by the field originating from the point 
source at {xo,yo) in the object plane. In the second in- 
terpretation, the first-zero full-width of the Airy pattern 
Co — 1.22Xzo/a represents the area of the object plane 
that most contributes to the field at the observation point 
{xi,yi) in the image plane. In the non-isoplanatic situa- 
tion described by Eq. (j4|), important interferences occur 
when the phase factor exp [i^r] cannot be considered as 
constant over the circular region of area 7r(eo/2)^ centered 
on {xi/M,yi/M) in the object plane. 

It is interesting to note that isoplanatism is recovered 
when one considers the imaging problem from So to Si (see 
Fig- HI); where So (Si) is a spherical surface of radius Zo 
(zi) centered on the pinhole, and having its vertex on O 
(X). The spherical field distributions on So and Si are de- 
scribed by the amplitudes Us„{xo,yo) and Usi{xi,yi), re- 
spectively (z-coordinates are dependent variables). These 
are related to the field amplitudes on O and X by a phase 
relation: Us,{xo,yo) = t/o(a;o, J/o) exp[i27r(r - Zo)/A] and 
UsA^i^Vi) = f/i(a;i,yi)exp[-i27r(s-Zi)/A]. Inserting these 
relations into Eq. ^ and removing the constant phase fac- 
tor exp[i27r(zo + Zi)/\] one gets: 
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which is a convolution relation, as expected for an iso- 
planatic imaging system. The physical reason why isopla- 
natism is recovered when the object lies on the spheri- 
cal surface So is simple to understand: since all the point 
sources are at same distance Zo from the pinhole, their point 
spread functions always interfere constructively. Note that 
there is no real need to measure the image on Si for it only 



differs from the image in the plane Z by a space-dependant 
phase that any standard intensity detector is insensitive to. 

3. Non-isoplanatism of a thin lens 

The scenario exhibited in the simple example of the pin- 
hole camera actually occurs in all optical imaging systems. 
When imaging with lenses or mirrors, however, the diffrac- 
tion effects are less dramatic than with a pinhole camera. 
The problem of phase distortion remains though, and may 
sometimes induce unwanted intensity modulation in the 
image. 




Fig. 3. Scheme of coherent image formation by an aberration-free 
thin lens. The points P^ and P^ in the image plane X are the geo- 
metrical images of the point objects Pi and P2 in the object plane 
O. The light-gray zones around P[ and Pj i^i the X-plane represent 
the individual diffraction patterns resulting from illumination of the 
pinhole by Pi and P2 respectively. The dark-gray zone represents the 
area where the individual diffraction patterns overlap and interfere. 

An aberration-free thin lens is an ideal diffraction-limited 
optical system that acts locally on the impinging field as a 
pure phase transparency r (a;, y) = exp[— z7r(a::^-|-2/^)/(A/)], 
where / is the focal length of the lens. The impulse re- 
sponse function of a thin lens |4|2j is an important concept 
that has proven to be very useful for the design of optical 
systems and optical data processing. It turns out that it 
is given by the same formula as for the pinhole camera — 
Eq. ([3]) — with the understanding that Zi is now related 
to Zo through the lens law l/z^ -I- l/z^ = 1//. To stress 
the analogy between the pinhole camera and the thin lens 
systems, we have drawn on Fig. [3] the same information as 
on Fig. [TJ and used the same notation. In particular, note 
that 2a now represents the diameter of the lens, which is 
usually considerably larger than the aperture of the pin- 
hole camera. As a consequence, the diffraction effects will 
be weaker, and non-isoplanatism will be less pronounced. 
Apart from this comment, the results and discussion of the 
last section also apply to the thin lens system. The relation 
between the field distribution in the plane O and its image 
on the plane T is non-isoplanatic and is given by Eq. Q . As 
with the pinhole camera, isoplanatism can be recovered — 
Eq. (O — when the thin lens is used to image the spherical 
surface So to Si. However, this is much less obvious here 
because, in contrast with the pinhole camera, a lens does 
not have an infinite field of view. That Si is the image sur- 



face of So is therefore questionable. This point requires a 
more carefull argumentation that we postpone to Sec. [5] 

Whether non-isoplanatism leads to interference when a 
plane emitter is imaged to a plane receptor depends not 
only on the size of the lens, but also on the size of the 
object and its position with respect to the optical axis. 
The condition for avoiding any interferences due to non- 
isoplanatism is that the phase 27rr/A in Eq. (|4]) varies by less 
than TT when the point corresponding to coordinates {xo , Ho) 
explores the Airy pattern (5_2_ [xi/M — Xo, yi/M — yo) in the 
object plane. If this is the case, the phase factor exp[i27rr/A] 
can be taken out of the integral in Eq. ((4|) . Since the fastest 
phase variation occurs when the point moves radially off 
axis, the following condition is obtained: 



.o(pg + ^ 



< XZo 



(6) 



where p% = ^J{xi/MY -f {yi/MY is the off-axis distance of 
the geometrical object point under consideration and Co = 
1.22Azo/a is the first-zero full width of the Airy pattern 
in the object plane. Note that Tichenor and Goodman [4] 
found a similar condition following a different reasoning. 
The main conclusions are: 

- When imaging points very close to the optical axis (pg <S 
60/4), interference due to non-isoplanatism does not oc- 
cur if a? /Xzo ^ 1. This last condition is always satisfied 
in practice with lenses. (For a diffraction-limited pinhole 
camera it fails to be satisfied.) 

- When imaging off-axis points (p^ ^ eo/4), interference 
due to non-isoplanatism does not occur if p^ ^ a, i.e. if 
the object points are not as far off-axis as the edges of 
the lens. 

When criterion ^ is satisfied for any point on the object, 
Eq. ([4]) can be written as 
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(7) 



where the phase factor comes from the second order approx- 
imation of s(xi, yi) + r{xi/M, yi/M) and constant phases 
have been removed. An intensity detector in the plane T 
will record the same image as in the isoplanatic case — 
Eq. ([5]). It should however be noted that the impulse re- 
sponse function is still raon-isoplanatic because the phase 
curvature has not been removed. It is important to keep 
this in mind when a phase-sensitive detector (hologram) is 
used and/or if further optical processing is needed. 

4. Non-isoplanatism and large field imaging 

According to the discussion in Sec. [31 non-isoplanatism 
has no effect on the intensity detected in the image plane 
if condition © is satisfied. In that case, Eq. ([7]) can be 
used instead of Eq. (|4|). In practical applications of coherent 
imaging, the assumption is usually made ^2 that the field 



mapping from the object space to the image space is given 
by Eq. ([7]) for any lens of a given optical system. This as- 
sumption is very convenient from a theoretical point of view 
since it makes the analysis of optical systems much sim- 
pler by the use of standard Fourier optics methods. From 
a technical point of view, however, good quality lenses (as- 
pheric lenses, for instance) usually have a small diameter 
because of manufacturing constrains, and care must there- 
fore be taken when implementing optical systems based on 
Eq. ([7|). In this section, we show that the slight difference 
between Eqs. ^ and ^ may lead to strong effects if the 
object does not fulfill condition ^. 

Let's first point out what is wrong with Eq. ^ from a 
physical point of view. Consider an object field Uo{xo, Ho) 
that is varying slowly on the length scale Co = 1.22Xzo/a of 
the peaked function 5a/(\z„){xi/M — Xo,yi/M — yo). Then, 
according to Eq. ([7]), one can write 
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This equation implies that the energy is conserved: 
\\\Ui{xi,yi)\^ Axidyi = JJ |C/o(a:o,2/o)|^ da;,, dj/o. However, 
independently of how slowly the field varies in space, 
Eq. ([8]) cannot hold for far off-axis points when the lens has 
a limited aperture. Radiation from far off-axis points (like 
P2 on Fig. [3|) is partially lost, and energy cannot be con- 
served. For instance, in the case of a plane wave travelling 
along the optical axis, radiation from the neighborhood of 
P2 will not be transmitted at all. We can therefore con- 
clude that Eq. ([7]) does not properly account for energy 
loss due to the limited aperture of the lens, especially for 
light originating from far off-axis points. Taking the phase 
factor exp(z27rr/A) out of the integral sign in Eq. ^ breaks 
down the energy balance. 

To give a deeper insight of the effect of the phase factor 
exp(i27rr/A), we use Eq. ([4]) instead of Eq. ^ to compute 
the intensity distribution in the image plane for a slowly 
varying object field. We now obtain 
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instead of Eq. ([8]). The paraxial approximation r ^ zq + 
{xl + yl)/{2zo) has been used. Making a change of integra- 
tion variables from (xq, j/o) to (C, 77) = (a;^, yo) - {x^,yi)/M, 
the preceding equation can be written as 
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The integral in Eq. ([5]) has the form of a Fresnel diffraction 
integral. For points {xi,yi) which are sufficiently far away 
from the optical axis, the quadratic phase terms C^/2 and 
77^/2 can, as a first order approximation, be neglected, and 
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Fig. 4. Intensity in the image plane X when a large square aperture, 
centered on the optical axis, is imaged with a thin lens. The square 
aperture is illuminated by a plane wave travelling along the optical 
axis. The sides of the square are b = 9.5 mm. The different panels 
correspond to different lens diameters. Prom left to right and top to 
bottom, the lens diameter is progressively increased in 2 mm-steps 
from 2a = 6 mm to 2a = 22 mm. The relevant parameters are: 
A = 780 nm, / = 12 mm, and Zo = 2 m. The scale of the figures is 
expressed in microns. The intensity scale is such that 1 corresponds 
to the expected uniform intensity in the center of the square in the 
limit of an infinite lens. 

the integral reduces to the Fourier transform of the Airy 
pattern in the object space. Using 

JJg,2.(.C+,,) ^^(^^^) ^^^^ ^ ^-^^ ^^^ j_\^ (10) 

Eq. ^ becomes 
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In Eqs. (fTO l) and pT | . 



circ(x, y) 



1 if ^/x^ + y^ < 1, 
if \/a;2 + j/2 > 1. 



(11) 



(12) 



Eq. (fTTjl shows that the intensity in the image plane exhibits 
a cut-off. No intensity reaches the image plane at a distance 
higher than Ma from the optical axis. This can be under- 
stood in the following way: since the field has been assumed 
to be slowly varying, the diffraction in the propagation from 
the object plane to the lens is negligible and the limited 
aperture of the lens has the same effect as a stop of radius a 
in the object plane. The function circ {xi/{Ma), yi/{Ma)) 
is the image of that virtual stop and can be interpreted as 
the shadow of the lens. This is, however, only a first order 



approximation, since the quadratic phase terms C^/2 and 
ry^/2 in Eq. ([9]) have been neglected. The effect of these 
quadratic phase terms is to create radial intensity oscilla- 
tions in the image, especially around the cut-off radius Ma. 

In Fig. m the preceding discussion is illustrated with an 
example. A square object of size 9.5 x 9.5 mm^ is imaged 
using a lens of diameter 2a varying from 6 mm to 22 mm. 
The object can be seen as a plane screen with a square aper- 
ture in it which is illuminated by a plane wave propagating 
along the optical axis. The object field Uo{xo, Vo) is 1 inside 
the square and zero outside of it. Fig.[4]shows the intensity 
distribution in the image plane computed using Eq. (|¥|). 
With small lenses (2a up to 10 mm) the clipping predicted 
by Eq. (fTTjl is observed. One can clearly distinguish the disk 
(fT2|) that limits the observable part of the object, as well 
as the intensity ripples due to the quadratic phase terms in 
© which were neglected when deriving pT|) . Interestingly, 
the intensity oscillations do not disappear as soon as the 
lens becomes bigger than the object. For 2a — 22 mm, some 
residual modulation still remains. Note that the parame- 
ters of this simulation are realistic ones: the object distance 
Zo has been fixed to 2 m (close to infinite-conjugate ratio 
imaging) and the focal length / to 12 mm (the numerical 
aperture ranges from 0.24 to 0.67). Note that energy is lost 
when imaging large objects through small lenses: For the 
simulations shown in Fig. [H the percentage of transmit- 
ted energy is, from left to right and top to bottom, 31.1%, 
56.1%, 83.4%, 96.6%, 99.2%, 99.8%, and nearly 100% for 
the last three images. 

The discussion leading to Eqs. ([8]) to pT|) only concerned 
object fields Uo{xo, Ho) that are slowly varying on the length 
scale Co — 1.22 Xzo/a in the object plane. For quickly vary- 
ing fields the previous discussion does not hold, but non- 
isoplanatism still has some effects on imaging. The way 
non-isoplanatism modifies the intensity distribution in the 
image plane strongly depends on the object wavefront. No 
general features can be drawn in that case. To get some 
insight, consider the following example. 

A large square grid (19 x 19 mm^) of mutually coherent 
point sources is imaged through a lens of diameter 2a = 
8 mm. The imaging conditions are otherwise the same as in 
Fig-Hi Let's consider that all the point sources are in phase 
and discuss the image formation when the spacing d be- 
tween the point sources is varied. If the points are well sep- 
arated {d ^ eo/2, so that the Airy patterns associated with 
them in the image plane do not overlap) no interference 
takes place and non-isoplanatism has no effect on the in- 
tensity distribution in the image plane. If the Airy patterns 
overlap, interferences similar to those depicted in Fig.[2]oc- 
cur. Two cases must be distinguished: d <^ eo/2 and d « 
eo/2. For d <§; eo/2, the object field varies slowly in space; 
a fringe pattern similar to the one in Fig. |4] is expected. 
The interesting case is d w eo/2. Fig. [5] shows the intensity 
distribution in the image plane when d ranges from 0.6 to 
2.2 times eo/2. For d < eo/2, a fringe pattern similar to 
the one in Fig.|4]is seen in the center of the field. However, 
periodic replicas of this pattern are also observed. For d > 
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Fig. 5. Intensity in the image plane I when a large squared 
(19 X 19 mni'^) grid of point sources, centered on the optical axis, 
is imaged by a thin lens. The point sources are mutually coherent 
and in phase. The different panels correspond to different values of 
the grid period d. From left to right and top to bottom, d is var- 
ied from 0.6 eo/2 to 2.2 eo/2 {eo/2 is the separation corresponding 
to the Rayleigh resolution criterion). The relevant parameters are: 
a = 4 mm, A = 780 nm, / = 12 mm, and Zo = 2 m. The scale of 
the figures is expressed in microns. The intensity scale is arbitrary 
but the relative dynamical range is indicated by the number in the 
upper right corner of each panel. For instance, the number "6" in 
the first panel means that the intensity range is 6 times larger than 
in the three panels of the last row. 

eo/2, the circular fringe patterns intersect each other, but 
the sparse sampling due to the grid structure of the image 
makes this structure barely visible (periodicity however re- 
mains). For clarity, only the central 30 x 30 pm^ region of 
the image is displayed in the last six panels of Fig. \5[ The 
simulations of Fig. [5] show that the intensity distribution 
in the image plane exhibit two distinct pseudo-periods (in 
both X and y directions) : the small-scale pseudo-period Md 
due to the grid structure of the object and the large-scale 
pseudo-period X associated with the ring patterns due to 
non-isoplanatism. Strictly speaking the image is periodic 
only if X is an integer multiple of Aid, in which case the 
period of the image is X. Hereafter, we use this property 
to deduce the value of X. The object field is modelled as a 
two-dimensionnal Dirac comb: 

Uo{xo, Vo) = 22 ^^^° - n d,yo-m d). 



where the integers n and m run from — oo to +cxd. Using 
Eq. ([U, we then find that 
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In order to find X, we require that \Ui{xi,yi)\'^ = |C/i(a;i + 
X, yi)p whenX is a multiple of Md. Using eq. (fT3|) . we have 
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The argument (xi+X) /M — nd, can be rewritten as Xi /M — 
n'd with n' ^ n — X/{Md) G Z. Replacing the sum over n 
by a sum over the values of n' , we obtain 
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One can note that 



2n' 
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Md 



and 



t>nm just differ by a factor 
n' X 2tt, and therefore \Ui{xi,yt)\'^ = \U^{xi +X,yi)\'^, if X 
is given by 



X = M X 



ZnX 



(14) 



When X is not a multiple of Af d, formula (|14|1 is still valid 
if X is understood as the large-scale pseudo-period. How- 
ever, the notion of "large-scale pseudo-periode" only holds 
if X/{Md) w a eo/d^ > 1 (e^ = 1.22 \Zo/a). Since a/d 
is usually a large number, this condition can still be valid 
even if d strongly exceeds the distance eo/2 corresponding 
to the Rayleigh criterion. The periodic patterns displayed 
in Fig.[5]prove that the field in the image plane correspond- 
ing to the Airy pattern of a given point in the grid influ- 
ences the entire image on a length scale much longer than 
the usually considered "Airy pattern diameter" Co- Readers 
familiar with signal processing will notice that the reason 
why periodicity appears here is the same one that makes 
the Fourier spectrum of a sampled signal (Dirac comb) pe- 
riodic. We however stress that the periodicity that is de- 
scribed here appears in the image plane and not in the 
Fourier plane of the object field. This is a peculiarity of 
non-isoplanatic imaging. 

Though non-isoplanatism can strongly influence the im- 
age formation of both slowly varying and rapidly vary- 
ing object fields, the effects are qualitatively different. For 
slowly varying fields, we have seen that the major effect 
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Fig. 6. Intensity in the image plane X when squared and circular 
contours are imaged by a thin lens. In the object plane, all the point 
sources forming the contours are mutually coherent and in phase. 
The imaging conditions are the same as in Fig. |4] (A = 780 nm, 
/ = 12 mm, and Zo = 2 m), but the size of the lens is fixed: 
2a = 6 mm, as in the first pannel of Fig. |4] The intensity scale is 
arbitrary. Panel (a) shows the image of three squared contours of 
increasing side b = a, b = 2a and b = 9.5 mm fa 3.2 a. The size of 
the outer square is the same as in Fig. |4] Panel (b) shows the image 
of three concentric circles, the diameters of which are equal to the 
sides of the squares in the panel (a). 



comes from the shadow that the lens projects in the image 
plane. The image is clipped so that only the central disk of 
radius Ma remains (see Fig. |4]). The simulations of Fig. \5\ 
show that, for rapidely varying object fields, light is ob- 
served outside the central disk of radius Ma. More complex 
effects are exhibited, but clipping do not occur anymore. In 
many applications of coherent optics (like optical lithogra- 
phy), the object is made of lines instead of single points or 
filled surfaces. Lines are objects on which the field varies 
slowly in one direction (the direction tangeant to the line) 
and rapidly in the orthogonal one. Therefore, one can ex- 
pect that the effects of non-isoplanatism will be intermedi- 
ate between the two previous cases. Fig. [6] helps to under- 
stand how lines are imaged through a non-isolanatic optical 
system. The imaging conditions are the same as in Fig. HI 
except that the lens diameter is fixed: 2a = 6 mm. Panel 
(a) shows the image of three squared contours. The outer 
square is exactly the contour of the filled square imaged in 
Fig-Hl The comparison with the upper-left pannel of Fig. [5] 
shows that the major part of the contour is now visible; 
only the corners of the square are clipped. Closer examina- 
tion shows that, for a straight line, only a segment of length 
2a in the object plane (2Afa in the image plane) is trans- 
mitted. The part of the Hne that is clipped corresponds to 
the orthogonal projection of the shadow of the lens on the 
straight line. This can be easily understood by analysing 
the propagation of the cylindrical waves emitted by straight 
lines through the spherical lens. The smaller square in pan- 
nel (a) is transmitted because its side is shorter than the di- 
ameter of the lens. The intermediate square is at the limit of 
the cut-off. As shown in the pannel (b) of Fig. [SI circles are 
never clipped, whatever their radii, because the orthogonal 
projection of the lens disk on the cicle is the circle itself. 



5. Isoplanatic imaging through a thin lens 

In Sec. [31 we claimed that isoplanatic imaging is possible 
with a thin lens when the object lies on the spherical surface 
So and the image is observed on the spherical surface Si. 
Let's examine this statement more closely. 

Referring to Fig. [3l simple Gaussian Optics arguments 
suggest that if the point P2 is translated horizontally from 
the object plane O to the surface So, its image P2, instead 
of moving towards Si , should move away from the lens. This 
argument is correct, but it relies on the Gaussian approxi- 
mation that the point P2 is initially in the "image plane" 
X. In reality, because of the field curvature aberration due 
to the lens, the stigmatic image of P2 is closer to the lens 
than the surface Si itself (Pj lies on the so-called Petzval 
surface) . Bringing P2 on So will place its image P2 exactly 
on Si [11]. Consequently, Eq. ([5|) is exact, while Eq. Q 
is only valid in the context of Gaussian approximation. If 
the lens is still diffraction-limited in this regime (no point- 
aberrations), the isoplanatic imaging geometry provides a 
nearly perfect transfer of the coherent field from the object 
to the image space. Only the resolution is reduced because 
of the finite size of the lens. 

The problem that remains is how to transfer a given 
field from a plane surface to a curved one before imaging, 
and vice versa after imaging. Generally, the object field 
is obtained by modulating a coherent illumination field 
A{xo, Ho) by a complex amplitude function Uo{xo, Ho)'- 



UoiXo^Vo) = A{Xo,yo) X UoiXo^Vo)- 



(15) 



The modulation Uo{xo, Ho) can be produced by any plane 
modulation device like a plane transparency, a grating 
or some kind of adaptative optics. Usually Uo{Xo,yo) it- 
self is considered as the "object" of interest. In that case, 
the illumination beam must be a plane wave propagating 
along the optical axis in order to map the profile Uo{xo, Ho) 
to the beam: Uo{Xo,yo) oc Uo{xo,yo)- If, instead, we il- 
luminate the modulation device with a spherical wave 
A{xo,yo) = Aexp (— z27rr/A) focusing on O (see Fig. [3]), 
we get UsA^o,yo) = ?7o(a;o,yo)exp(j27rr/A) oc Uo{xo,Vo)- 
This kind of illumination projects the object onto So as 
required for isoplanatic imaging i2]. Starting with an illu- 
minating plane wave travelling along the optical axis, this 
can be achieved by placing an additional thin lens of focal 
length Zo just before (or just after) the object plane O. 
This lens must act as a pure phase-correction transparency. 
Similarly, on the image side, a thin lens of focal length Zi 
placed just after (or before) the image plane can be used 
to project the image from Si onto I. From a broader point 
of view, any spherical field distribution in the object space 
of a centered paraxial optical system can be imaged onto a 
spherical surface of any curvature using only lenses; some 
of them will be imaging lenses, while others will play the 
role of phase-correction transparencies. This is the basis 
of the so-called metaxial optics theory formulated by the 
Bonnet |12|13|14j . 



It should be noted that this approach only works for suf- 
ficiently slowly varying fields, because the diffraction from 
O to So (and Si to X) has to be negligible for the amplitude 
of the fields on O and So {Si and 2) being the same. In ad- 
dition, due to optical design constraints, the use of optical 
lenses as phase-correction transparencies is not always pos- 
sible. Moreover, lenses are never perfectly thin and diffrac- 
tion through them may have worse effects on imaging than 
the phase curvature due to non-isoplanatism. Whether or 
not non-isoplanatism should be corrected depends on the 
particular system under consideration. A clever design can 
minimize, if not cancel, its effects on the detected intensity. 



6. Experimental investigation 

The theoretical discussion of Sees. |4] and [5] relies on two 
strong approximations: the paraxial approximation and the 
thin lens approximation. One may wonder whether our 
analysis is robust enough to be applied to systems contain- 
ing powerful lenses, which are usually thick and have a high 
numerical aperture. The following experiment shows that 
the previous discussion is also valid for these systems. 

The setup is shown in Fig. \7\ The test lens Li is an as- 
pheric lens having a focal length / = 8 mm and a diam- 
eter 2a = 8 mm (LightPath 352240). This lens is used to 
image the surface of a spatial light modulator (SLM), a 
1024 X 768 micromirror array, that modulates the ampli- 
tude of the reflected beam. Using the SLM, we can gener- 
ate arbitrary patterns. The resolution is set by the size of 
the micromirrors (13 x 13 pm^). The object plane is 60 cm 
away from the lens. We use a double- lens system (L3, L4) 
to magnify the image produced and project it on the CCD 
camera. L3 is a diffraction-limited aspheric lens and L4 a 
long- focal achromatic doublet. The numerical aperture of 
this double- lens system is large enough to prevent any pos- 
sible clipping or diffraction during the magnification pro- 
cess. The SLM is either illuminated with a plane wave or a 
spherical wave converging on Li (A = 780 nm). In the first 
case, the imaging system is non-isoplanatic. As explained 
in Sec. [5l it becomes isoplanatic when a spherical-wave il- 
lumination is used. The spherical wave is obtained from 
the impinging plane wave by inserting the additional lens 
L2 (75-cm focal length, achromatic doublet) in front of the 
SLM. We use this setup to image objects of different sizes 
and shapes through the test lens. 




Fig. 7. Experimental setup: Li, aspheric lens (/ = 8 mm); L2, achro- 
matic doublet (/ = 750 mm); L3, aspheric lens {/ = 20 mm); and 
L4, achromatic doublet (/ = 500 mm); SLM, spatial light modula- 
tor; CCD, coupled-charge camera; O, object plane; I, image plane. 





Fig. 8. Images of different-sized squares (the side b ranges from 
1.4 mm to 9.6 mm) recorded using the setup in Fig. [7] The object 
plane is either illuminated with a parallel beam (non-isoplanatic case, 
left column) or with a spherical one converging on Li (isoplanatic 
case, right column). 



Fig. 9. Images of a rectangular grid of points recorded using the 
setup in Fig. [7] The object plane is illuminated with a parallel 
beam (non-isoplanatic case). The period of the grid in the object 
plane, d = 109 jam, corresponds to 1.5 times the Rayleigh criterion 
separation (eo/2). 

Fig. [8] shows the pictures recorded by the CCD camera 
when we image squares of different sizes through the 8-mm 
diameter aspheric lens Li . Let's first consider the case when 
the SLM is illuminated by a parallel beam (non-isoplanatic 
imaging, left column of Fig. [8]). For a square side larger 
than 8 mm, the clipping effect described in Sec.|4]is clearly 
observed. Circular fringes similar to those of Fig.[4]are also 
seen in the image plane. For a square side smaller that 
8 mm, the region of non-zero intensity is limited by the size 
of the square. However, intensity modulation due to non- 
isoplanatism is still noticeable for a square side as small 
as 2.7 mm. In the case of the 1.4 x 1.4 mm^ square, non- 
isoplanatism has a negligible effect. Note that, in the sim- 
ulations of Fig.m the lens size was varied while the object 
size was kept constant. Here, the lens is always the same, 
but the square size is varied instead. For this reason, the 
ring pattern is the same for every picture in the left column 
of Fig. [H When the SLM is illuminated with a spherical 
wave converging on Li (isoplanatic imaging, right colomn 
of Fig. [8]) no clipping effect occurs and there is, in principle, 
no limit to the size of the objects that the system can im- 
age (the slight variations in intensity are due to amplitude 
inhomogeneities is the illumination beam). 

In a second experiment, we imaged rectangular grids of 
points illuminated by a plane wave (non-isoplanatic illu- 
mination), a situation that we described theoretically in 
Sec. m Fig. [9] shows the recorded intensity distribution in 
the image plane for a grid with a period d — 109 ym, which 
corresponds to 1.5 times the Rayleigh criterion separation 
(eo/2). As shown in Sec. IH for such large value of d the 
ring patterns should not be visible anymore (see the six 
last panels of Fig.[5|). Here, however, we can distinguish one 



central ring pattern and four replicas intersecting it. We 
attribute this discrepancy to the fact that the experimen- 
tal point-spread function is a bit broader than the theoret- 
ical Airy pattern. As noted is Sec. IH the "disappearance" 
of the multiple ring patterns when d increases is caused by 
the sparseness of sampling due to the grid. A broadening 
of the point-spread function reduces that sparseness and 
therefore restores the ring patterns characteristic of non- 
isoplanatic imaging. 

7. Conclusion 

In general, even diffraction-limited imaging systems dis- 
tort the phase of the processed fields. This is of no relevence 
when used with incoherent light, but has a tremendous ef- 
fect on coherent imaging. In combination with Fraunhofer 
diffraction from the finite instrumental aperture, the phase 
distortion leads to a severe degradation of the field ampli- 
tude in the image plane. We analyzed this phenomenon for 
two very different but prototypic imaging systems (the pin- 
hole camera and a thin lens) and observed its fundamental 
and general nature. We showed that substantially different 
effects arise depending on whether the field varies slowly or 
rapidly on the length scale of an Airy pattern. The degra- 
dation of the field amplitude can however be overcome or, 
at least, limited by a clever design of the optical system. 
The main aspects of our analysis have been confirmed ex- 
perimentally using a powerful thick aspheric lens to demon- 
strate that the phenomenology that is described also holds 
beyond the paraxial and thin lens approximations used for 
the theoretical analysis. 
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